#include "Cphase.h"

Cphase::Cphase()
{
	V = vector<double>(2, 0.0); P = 0.0; rho = 0.0;
	T = 0.0; mu = 0.0; kappa = 0.0; cp = 0.0; gamma = 0.0; Pr = 0.0; Sc = 0.0; Dva = 0.0;
}

Cphase::Cphase(const Cphase& c0)
{
	V = c0.V; P = c0.P; rho = c0.rho; T = c0.T; mu = c0.mu; kappa = c0.kappa;
	cp = c0.cp; gamma = c0.gamma; Pr = c0.Pr; Sc = c0.Sc; Dva = c0.Dva;
}

Cphase::Cphase(const shared_ptr<CaseSettings> caseSetting1, const shared_ptr<Constants> const1)
{
	V = vector<double>(2, 0.0);
	P = (*caseSetting1).inletP;
	T = (*caseSetting1).inletT;
	mu = (*caseSetting1).inletMU;
	rho = (*caseSetting1).inletRho;

	kappa = (*const1).k_air;
	cp = (*const1).Cp_dryAir;
	gamma = (*const1).gammaAir;
	Pr = mu * cp / kappa;
	Sc = mu / (rho * (*const1).D0);
	Dva = (*const1).D0;
}